Isospin-symmetry restoration within the nuclear density functional theory: formalism 

and applications 
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Isospin symmetry of atomic nuclei is explicitly broken by the charge-dependent interactions, pri- 
marily the Coulomb force. Within the nuclear density functional theory, isospin is also broken 
spontaneously. We propose a projection scheme rooted in a mean field theory, that allows the con- 
sistent treatment of isospin breaking in both ground and exited nuclear states. We demonstrate 
that this scheme is essentially free from spurious divergences plaguing particle-number and angular- 
momentum restoration approaches. Applications of the new technique include excited high-spin 
states in medium-mass N—Z nuclei, such as superdeformed bands and many-particle-many-hole 
terminating states. 



PACS numbers: 21.10.Hw, 21.60.Jz, 21.30.Fe, 27.40. 



I. INTRODUCTION 



Density functional theory (DFT) in its original formu- 
lation [f|, 0] and its extended versions 0, 0, Ha, 0] has be- 
came a universal approach to compute the ground-state 
(g.s.) and excited [7], Q configurations of many-electron 
systems held together by an external one-body potential 
in condensed-matter, atomic, and molecular physics. The 
DFT strategy has also been used in the area of nuclear 
structure. The nuclear DFT, a natural extension of the 
self-consistent mean- field (MF) theory d,[T| is a tool of 
choice for computations of g.s. properties and low-lying 
excitations of medium-mass and heavy nuclei. 

There are a number of differences between the elec- 
tronic DFT and nuclear DFT. First of all, because 
of the absence of external potential, nuclei are self- 
bound systems, and this creates conceptual problems due 
to the distinction between intrinsic (symmetry-broken) 
and laboratory-system densities. This problem can be 
treated within the extended DFT framework of Levy- 
Lieb 0, 1, 1, 1 p, see Refs. [H 0, 0, El. The 
second difference is related to the effective interaction. 
A complicated form of the effective nucleon-nucleon in- 
teraction, as well as complexity of the nuclear matter 
saturation- mechanism (pertaining to three-body forces), 
give rise to the presence of higher-order terms in the nu- 
clear energy density functional (EDF). Hence, the EDF 
cannot be constructed in a model-independent ab initio 
way based solely on bulk properties of infinite-medium 
augmented by exchange and gradient corrections, as it 
is usually done for electronic systems. The current best 
nuclear EDFs are still phcnomcnological in nature, ei- 
ther inspired by a MF averaging of the effective interac- 
tion such as the zero-range Skyrme i6| force, or based 
on systematic expansion involving symmetry-constrained 
gradient terms [TtJ- By construction, the nuclear EDF 
is characterized by a certain number of free parameters 
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which are fitted directly to empirical data. 

The third major difference between electronic and nu- 
cleonic DFT, addressed in this study, is the presence of 
two types of fermions and the charge independence of 
the nuclear interactions giving rise to the isospin symme- 
try. Since the isospin symmetry is violated essentially by 
the Coulomb interaction, which is much weaker than the 
strong interaction between nucleons, many effects asso- 
ciated with the isospin breaking in nuclei can be treated 
in a perturbative way, making the formalism of isotopic 
spin a very powerful concept in nuclear structure and 
reactions 0, [f9j]. 

Up to electromagnetic effects, the isospin symmetry 
should be conserved by elementary excitations of nuclei. 
This is not the case for the elementary excitations modes 
of the nuclear Hartree-Fock (HF) or Kohn-Sham theory, 
i.e., proton and/or neutron particle-hole (p-h) excita- 
tions. The independent-particle wave function manifestly 
breaks the isospin symmetry in the g.s. configurations of 
odd-odd N—Z nuclei and in all other but isoscalar ex- 
cited states. Two prominent examples are discussed in 
this work: superdeformed (SD) bands in a doubly magic 
nucleus 56 Ni [20| and terminating states in N=Z A ~40 
nuclei [Hj]. 

The paper is organized as follows. We begin in Sec. |TT] 
with general discussion of g.s. configurations of odd-odd 
N=Z nuclei and p-h excitations in even-even (c-c) nu- 
clei, so as to qualitatively introduce problems faced by 
MF-based theories around N=Z and to provide motiva- 
tion for the necessity of isospin restoration. Section Ifffl 
contains a detailed presentation of the isospin-projected 
DFT approach. In particular, we demonstrate that this 
approach is essentially free from problems related to un- 
compensated poles plaguing projection techniques [22| : 
hence, it requires no further regularization j23|. Applica- 
tions of the isospin-projected DFT approach to the SD 
bands in 56 Ni and terminating states in N—Z isotopes 
of Ca, Sc, Ti, V, and Cr, are presented in Sees. IfV Al 
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FIG. 1: The upper panel schematically shows two possible 
g.s. configurations of an odd-odd N=Z nucleus, as described 
by the conventional deformed MF theory. These degener- 
ate configurations are called aligned (left) and anti-aligned 
(right), depending on what levels are occupied by the valence 
particles. The lower panel shows what happens if the isospin- 
symmetry is restored. The aligned configuration is isoscalar; 
hence, it is insensitive to the isospin projection. The anti- 
aligned configuration represents a mixture of the T=0 and 
T=l states. The isospin projection removes the degeneracy 
by shifting the T=0 component down. 



and IIV Bl respectively. Finally, Sec. fVl contains the con- 
clusions of this work. 



II. ISOSPIN SYMMETRY VIOLATION IN A 
MEAN-FIELD APPROACH 

Let us begin by exposing serious problems with the 
MF description of g.s. configurations of odd-odd (o-o) 
N=Z nuclei. If, for the sake of simplicity, the Coulomb 
and time-odd polarization effects are disregarded and 
proton-neutron symmetry is conserved, a deformed MF 
approach naturally leads to four-fold degenerate (isospin 
and Kramers) single-particle (s.p.) levels. Consequently, 
the MF g.s. configuration of the o-o N=Z nucleus is not 
uniquely defined and depends on occupation of specific 
levels. As shown in the upper panel of Fig. [TJ the va- 
lence proton and neutron can be arranged in two distinc- 
tively different ways. Indeed, configurations of the type 
shown on the left-hand side (aligned configurations), are 
symmetric in spin-space coordinates and, therefore, anti- 
symmetric (isoscalar) in isospin coordinates. The anti- 
aligned configurations depicted on the right-hand side, 
have a mixed symmetry in spin-space coordinates and, 
therefore, also in isospin coordinates. 

Because of their isoscalar character, aligned configura- 
tions are not affected by the isospin projection. On the 
other hand, the projection lifts the degeneracy of T=0 
and T—l components of aligned configurations as it is il- 
lustrated in the lower panel of Fig.[TJ Due to the repulsive 
character of the nuclear symmetry energy, the isovector 
(isoscalar) components of the anti-aligned configurations 
are shifted up (down) in energy. Hence, isospin restora- 
tion changes the structure of the ground states of o-o 
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FIG. 2: The upper panel schematically shows the s.p. configu- 
ration corresponding to the SD Band 1 in 56 Ni, which is based 
on a 4p4h excitation with respect to the doubly magic core. 
The middle panel depicts two possible configurations of Band 
2 involving a single proton (left) or neutron (right) promoted 
to the Ogg/2 shell. To facilitate the discussion, we also show 
the relevant asymptotic Nilsson quantum numbers associated 
with deformed orbitals involved. The lower panel illustrates 
the isospin splitting induced by the isospin-symmetry restora- 
tion. 



N=Z nuclei without affecting the ground states of e-e 
N=Z nuclei. In other words, it does affect the binding- 
energy staggering along the N=Z line. 

Another example that nicely illuminates problems with 
isospin encountered in MF approaches is the case of two 
SD bands observed in the doubly magic nucleus 56 Ni [2(ij |. 
Following Ref. [20], we label them as Band 1 and Band 2. 
Band 1 is interpreted as a four-particle four-hole configu- 
ration formed by promoting two protons and two neu- 
trons from the O/7/2 shell to O/5/2 shell. Within de- 
formed MF model, where it is more natural to use the 
notion of Nilsson orbitals, Band 1 is obtained by emp- 
tying the neutron and proton [303] 7/2 Nilsson extruder 
orbitals and occupying the prolate-driving [321] 1/2 lev- 
els. This interpretation, not involving the 0<?g/2 shell, is 
strongly supported by both state-of-the-art shell-model 
(SM) calculations [H 0] in the Ohu fp space and the 
self-consistent cranked-HF theory [2(|. Both approaches 
appear to reproduce satisfactorily the excitation energy 
and moment of inertia (Mol) of this band. The p-h con- 
figuration of this band is depicted schematically in the 
upper panel of Fig. [2] 

Band 2, on the other hand, cannot be reproduced 
within the fp SM space. This band has also a larger 
Mol as compared to the Band 1, and it is becoming yrast 
around spin / ~ 12ft. These two facts strongly suggest 
that its structure must involve at least one particle in the 
prolate-driving [440] 1/2 Nilsson orbital that originates 
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from the Ogg/2 shell. According to the calculations of 
Ref. (2(J, Band 2 involves one proton in 7r[440]l/2 state. 
Its configuration, which can be regarded as one proton 
p-h excitation with respect to the reference Band 1, is 
schematically shown in Fig. [2] (middle panel, left-hand 
side). This scenario has been supported by the full pf— 
and p/59/2-SM calculations [25]. 

The conventional MF interpretation of Band 2 is not 
fully supported by experiment. Indeed, the MF theory 
predicts existence of the second band that is built upon 
the single neutron occupying the v [440] 1/2 orbital (mid- 
dle panel of Fig. [2 right-hand side). This neutron band 
is predicted to be placed only slightly higher in energy 
than the proton band, and this is due to a tiny difference 
in the Coulomb energies between these two configura- 
tions with, otherwise, very similar properties. Therefore, 
it is extremely difficult to understand within MF the- 
ory why only one of these two bands is observed exper- 
imentally. Moreover, HF calculations predict that the 
ir [440] 1/2 band is too high in excitation energy and that 
it crosses Band 1 above spin / ~ 16ft, that is, well above 
the empirical crossing. In these discussion, energetic ar- 
guments essentially exclude configurations involving two 
or more particles in the Ogg/2 shell. 

Here again the problems with the MF approach can 
be traced back to the isospin-symmetry violation by the 
p-h excitations. Indeed, the MF approach treats proton 
\iXp-ti) and neutron \v p -h) P _n excitations as independent 
elementary excitations, and this manifestly breaks the 
isospin symmetry in N=Z nuclei. To make the elemen- 
tary excitation modes in N—Z nuclei consistent with the 
symmetry, one needs to symmetrize or anti-symmetrize 
the Slater determinants corresponding to the mirror p-h 
excitations \T = 0(1)) s» -j=(\~K p -h) ± |z/ p _h,)). Such wave 
functions go beyond the usual MF picture. Restoration 
of isospin lifts the MF degeneracy between proton and 
neutron excitations and shifts the isoscalar (isovector) 
configurations down (up) in energy, as illustrated in the 
lower panel of Fig. [5] Of course, it is irrelevant whether 
the T=0 and T=l components are projected from the 
proton or neutron state; the results should be identical 
up to the tiny polarization effects. 

At present, it is not at all clear how to include the 
isospin correlations directly into the functional. The 
most natural method of taking them it into account is 
the isospin projected DFT approach [H, HI HI, 
which is also used in this work. We note that such an 
approach is within the provisos of the DFT, whereby 
the isospin-projected energy is still a functional of the 
isospin-unprojected density, cf. Refs. [30l. l3lT 1321] . 

We finish the discussion in this Section by noticing that 
even the isospin-projected DFT approach cannot fully 
account for the structure of the g.s. of an 0-0 N=Z nu- 
cleus. Indeed, within the conventional DFT approach, 
the aligned and anti-aligned configurations differ due 
to different time-odd (TO) polarization effects in these 
two configurations. This polarization affects the posi- 
tion of the T = 1, T z = states with respect to the 



T = 1,T Z = ±1 states. While the former states are 
shifted due to TO polarization, the latter ones - the 1 = 
ground states of e-e N — Z = ±2 isotopes - are not influ- 
enced by TO effects due to time-reversal symmetry con- 
servation. Hence, the TO terms introduce a very specific 
source of the isospin symmetry violation in the intrinsic 
system. 

Moreover, as shown in the lower part of Fig. [T] the 
isospin projected DFT approach always yields ground 
states with T=0 in 0-0 N=Z nuclei. This is at vari- 
ance with empirical data; indeed, it is well known (see, 
e.g., Refs. [H, [H, HH) that, with two exceptions, the 
isospin of the g.s. changes from T=0 in light (A < 40) 
0-0 N=Z nuclei to T=l in heavier (A > 40) 0-0 N=Z 
nuclei. To be able to address this question in the nu- 
clear DFT, one has to consider simultaneous isospin and 
angular-momentum projections and also consider resid- 
ual proton-neutron interactions between valence particles 
(either perturbatively or by breaking the proton-neutron 
symmetry). This subject will be considered in forthcom- 
ing studies. The focus of this paper is on high-spin su- 
perdeformed and terminating states, which - due to large 
spin polarization and high seniority - are expected to be 
less affected by the effects mentioned above. 



III. THEORY 

A. Isospin projection formalism 

Within the MF approximation, the isospin symme- 
try breaking has two sources. The spontaneous break- 
ing of isospin associated with the MF approximation it- 
self [2^, [29|, [36| is, in this theory, intertwined with the 
explicit symmetry breaking due to the Coulomb interac- 
tion. The essence of our method [H, [29[ is to retain 
the explicit isospin mixing. It is achieved by rediag- 
onalizing the total effective nuclear Hamiltonian -ff/viv 
in good-isospin basis. The basis is generated by using 
the projection-after-variation technique, that is, by act- 
ing with the standard one-dimensional isospin-projection 
operator on the MF product state | $) : 



\TT, 



1 



\/Ntt z 
2T+ 1 



TT Z JO 



dfi sin /3 £r zT Jp) #(/3) I $),(!) 



where (3 denotes the Euler angle associated with the 
rotation operator R{(3) = e _I/3Ty about the y-axis in 
the isospace, dj, T (/?) is the Wigner function [33], and 
T z = (N—Z)/2 is the third component of the total isospin 
T. The wave function |$) is obtained self-consistently by 
solving the Skyrme-Hartree-Fock (SHF) equations. In 
the HF limit, the z-component T z is strictly conserved. 
The overlap Ntt z > or interchangeably the normalization 
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factor &t,Tz ! are given by 



2T+ 1 



dPsmp& Ti {p)tf{p), (2) 



where 



(3) 



stands for the overlap kernel. The good-isospin basis is 
spanned by the states \T,T Z ) that have tangible contri- 
butions to the MF state 



T>\T Z \ 



>T.T Z 



\T,T Z 



E 

T>\T X \ 



>T.T r 



1. 



> 10 



(4) 
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In practice, we retain states having |&t,t 2 
which sets the limit of T < \T Z \ + 5. 

The total nuclear Hamiltonian consists of the kinetic 
energy T, the Skyrme interaction V s and the Coulomb 
interaction V that breaks isospin: 



H 



NN 



= T + V + V 



+ V+ V& + V& + Vr 



10 



"20-1 



(5) 



where 
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are, respectively, isoscalar, covariant rank-1 (isovector), 
and covariant rank-2 axial (isotensor) spherical tensor 
components of the Coulomb interaction. They are con- 
structed by coupling spherical components of the one- 
body isospin operator: 



no 



n±i =^4= (T x ±iT v ), (9) 



where fj, i = x, y, z denote Pauli matrices and symbol o 
stands for the scalar product of isovectors. 

The Coulomb interaction is the only source of the 
isospin symmetry violation in our model. Charge sym- 
metry breaking components of the strong interaction and 
the isovector kinetic energy (which is quenched as com- 
pared to its isoscalar counterpart by a factor ~ 0.001 
due to very small mass difference AM between the neu- 
tron and the proton relative to mean nucleonic mass M) 
are not taken into account. Hence, the isoscalar part of 
the total Hamiltonian reads Hoo = T + V S + Vq^. Its ma- 
trix elements can be cast into one-dimensional integrals: 



(TT z \H oa \T'T : 
S TT ' 2T + 1 



5tt' 
N- 



TT- 



($|jW%t,I*> 



N- 



TT- 



d(3 sin/3d£ T (0) Woo 03), (10) 



involving the Hamiltonian kernel 

Woo 09) = ($|£oo-R(/3)|$> 



(11) 



Calculation of matrix elements of the isovector ([7]) 
and isotensor (jSj) components of the Coulomb interac- 
tion is slightly more complicated. It appears, however, 
that these matrix elements can also be reduced to one- 
dimensional integrals over the Euler angle /?: 



(T'T Z \V°\TT Z [ 



fi—X 

E 



(12) 



z p.— — A 



where 



(_-\\nn T ' T z (~< T ' T : 



2T+ 1 



TTzXO^TTz-iiXii 
ra\ /^IT>C 



(13) 



dp sin/? 4,-^os) mv&Rm*), 



and symbols C^j^xu stand for the Clebsch-Gordan coef- 
ficients [37( | . The non-axial (n^O) isovector and isoten- 
sor components of the Coulomb interaction appearing in 





Eq. (HH) are: 


(6) 




(7) 




(8) 


^2±2 ( r ij) 



If! 

4r, 



O r i±i 



'i±i 

1 e ' 1 (Ai) Ai) 
4r« V3 V Tl±lTl ° 



(14) 



- W *^ » , (15) 



T io T l±l 



4r« V 3 T 



(i) Af) 
1±1 T 1±1- 



(16) 



Equations (JTSJ) — (JTHJ) can be derived by using the standard 
transformation rules for a covariant spherical-tensor op- 
erator T\f+ of rank-A under the three-dimensional rotation 
by the Euler angles = (a, 0, 7), that is, 

R(Q)f x ,R(ny = J2 ^(^TV, (17) 

A*' 

where the SO (3) projection operator P^ M (f2) commutes 
with spherical tensors as [38[ 



' K. M . ' V» I'm, K, 

EC-D^^^w^^- (18) 



The cornerstone of the isospin projection scheme de- 
scribed above is a calculation of the Hamiltonian and 
norm kernels and subsequent one-dimensional integration 
over the Euler angle ft. The integrals are calculated nu- 
merically using the Gauss-Legandre quadrature, which is 
very well suited for this problem provided that the cal- 
culated kernels are non-singular. 
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B. Inverse of the overlap matrix 

A prerequisite for the calculation of the integrals in 
Eqs. ((2]), (fTU|) . and (pHI)) is the isospin rotation of the 
Slater determinant. To this end, one has to perform an 
independent rotation of s.p. neutron, 



'•p. 



- V. 



» 



and proton, 



(P) 



i = 1,2, 



i = N+l,. 



,N, 



(19) 



(20) 



HF wave functions. The use of the two-dimensional 
(spinor) representation of isospin constitutes a natural 
formalism for the isospin projection technique. In this 
representation, the upper (neutron, q — 1) and lower 
(proton, q = —1) components of the isospin-rotated 
states are given by: 



R(J3) w = <fi(p), 



(21) 



class of functionals (or effective density-dependent inter- 
actions) solely involving integer powers of local densities 
[23l ]. Unfortunately, almost all commonly used Skyrme 



and Gogny parameterizations, except for SIII 41|, in- 
volve fractional powers of the density. The appearance 
of singularities prevents us from using the local Slater ap- 
proximation for the Coulomb exchange. In the present 
work we treat it exactly using the method of the Gaussian 
decomposition of the Coulomb interaction, as described 
in Refs. E2,IH. 



Compared to the particle number or the angular mo- 
mentum projection schemes, isospin projection is a rel- 
atively simple procedure. In particular, the dependence 
of the inverse matrix 0~ 1 (/3) on the isorotation angle 
(3 can be determined analytically and this enables us to 
demonstrate that the isospin projection is free from ker- 
nel singularities. To this end, we write the overlap matrix 
(125T) in the form: 



o{ i) = ( °™I In 



smfOt cosf/z 



sinfO 

2 



(26) 





(n) a 
tp\ 'cos% 

(™) • 

<Pi Sm 2 

-^ (p) sin I 
^cosf 



To calculate the kernel ($|Q|$) of an arbitrary op- 
erator Q between two non-orthogonal Slater determi- 
nants, we apply the generalized Wick's theorem, see, e.g., 
Ref. [39( | . In particular, the norm kernel can be written 
in a compact form: 



($|$) = Det O, 
where the overlap matrix is: 



(24) 



(25) 



To calculate kernels appearing in the projection for- 
malism, one needs to invert the overlap matrix O^ 1 . This 
can cause serious problems due to the presence of singu- 
larities [22|, H3|. The regularization of kernel singulari- 
ties is a difficult problem [23|. Thus far, a regulariza- 
tion scheme has been worked out only for a very specific 



where In (Iz) stands for the N x N (Z x Z) unit matrix, 
while O = {(p( 

»)|p(p)) is the rectangular N x Z overlap 
matrix of the neutron and proton s.p. wave functions. 
Using the singular value decomposition (SVD) technique, 
the matrix O can be written as: 



o = wdvK 



(27) 



where W is an N x Z rectangular matrix having orthog- 
onal columns (W^W = Iz), D is a diagonal 2xZ real 
and non- negative matrix (D* = D), and V is a Z x Z 
quadratic unitary matrix (VW = VV^ = I z ). The SVD 
decomposition Q27[) further implies that: 



OO^ = WD 2 W\ 
O^O = VD 2 VK 



(28a) 
(28b) 



where OO^ and O^O are matrices of dimension N x N 
and Z x Z , respectively. 

The SVD decomposition allows us to analytically di- 
agonalize the overlap matrix 0((3). Without loss of gen- 
erality we assume that N > Z. Hence, the product 
0(—/3)0((3) can be written as: 



J 



O(-0)O(P) 



W 




cos 



2 a 



sin 2 %D 



.■2 0; 



D 2 



W 




(29) 
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where W of dimension N x N is a unitary matrix com- 
posed of columns W: 

W = (W, W) where WW^ = W^W = I N , (30) 

and W is the unitary complement. 

The D matrix, on the other hand, is the Z x Z matrix 
D completed to the dimension TV x TV by zeros: 

*-(?•). (3D 
I 



This expression shows explicitly that the inverse of the 
overlap matrix 0(/?)r. x can be singular only if [3 = n. 
Indeed, in the case of N=Z, the ith denominator becomes 
zero, 

^cos 2 ^I N + sin 2 ^D 2 ^j = 0, (33) 

only when (3 — n and the corresponding singular value 
vanishes: Di=0. In the case of N ^ Z, the only differ- 
ence is that the matrix D contains zeros by construction, 
see Eq. (|3"Tj) . Fortunately, for the isospin projection, the 
singularity at [3 = n is compensated by the Jacobian 
sin (3 d(3 . 

C. Isospin structure of the density matrix 

To calculate the transition density matrix and deter- 
mine its isotopic structure it is convenient to introduce 
auxiliary ket-states defined as: 

A 

q,(3) = &0*> 0) 0(/3)7/. (34) 

3=1 

where the space-spin coordinates are abbreviated as x = 
(r,cr). This allows for rewriting the transition density 
matrix into a diagonal form: 

A 

p(x,q;x',q') = ^ip*{x' ,q') 4>i(x,q,0), (35) 

i=l 

where the summation runs over a set of s.p. states, in 
full analogy with the HF particle density matrix. To 
alleviate the notation, in what follows we do not indicate 
the dependence of the transition density matrices on the 
isorotation angle (3. 



Since the first and third matrices on the right hand side 
of Eq. (|29p are unitary and the second matrix is diagonal, 
the inverse of the overlap matrix reads: 



(32) 



I 

An arbitrary matrix in the isospace can be decomposed 
in terms of the isospin Pauli matrices. In particular, the 
transition density matrix p(x, g; x\ q') can be expressed 
as: 

p(x,q;x',q') = ^p 00 {x,x')5 qq , 

+ 2 T la (q,q')p la (x,x / ), (36) 

Q!=0,±1 

where 

no{q, q') = (q\haW) = ?<W , (37) 

n±i(<7,<?') = (q\h±iW) = TV25 q , ±1 S q ,^x. (38) 

The isoscalar poo and isovector p\ a transition density ma- 
trices are calculated by contracting p: 

poo = Tr p, (39a) 
pio = Tr (rio/5), (39b) 
pi±i = -Tr (ti t i/5). (39c) 

The use of the two-dimensional spinor notation is 
convenient for analytical derivations. However, from 
the numerical-programming perspective it is better to 
present expressions in a one-dimensional form, because 
it allows us to exploit the separability of proton and 
neutron wave functions and densities - a common fea- 
ture of essentially all the existing HF codes including the 
HFODD solver used in this work. 

To write the isoscalar and isovector densities in the 
one-dimensional form we introduce two sets of neutron- 
like (i = 1, 2, . . . , N) wave functions: 

N 

~^\x) = Y,<pf(*)°W7i, (40) 
i=i 

A 

$\x) = Yl <p?H*)d(0)£, (41) 

j=N+l 



0(0)- 



W 
V 



,2 



D 2 



cos 2 ^In 



,2 



W 
V 



O(-0). 
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and two sets of proton- like (i = N + 1, N + 2, . . . .A) wave 
functions: 

~^\x) ee f>f>(x)6(/?)7.\ (42) 
ee £ <pf{x)6{P)£. (43) 

J=iV+l 

Again, the dependence of these states on the isorotation 
angle is not explicitly indicated. Using these wave func- 
tions, we define two neutron-like densities: 

N 

p(™> (*,*') = X><">V)# W) (*) (44) 

i=l 

iV 

p(f)(x,x') = E^V)^*). (45) 

i=l 

and two proton- like densities: 

#W) - E l^V)^*) (46) 

i=AT+l 

#W) = E ^V)^*)- (47) 
1=^+1 

In these definitions, subscripts indicate whether the den- 
sity is neutron- or proton-like and superscripts indicate 
whether the summation indices in Eqs. (|40[) - (|43[) run over 

neutron or proton states. The densities p„\x,x') and 
pp(x,x') can be associated, respectively, with the rais- 
ing and lowering components of the analog spin [H, [45[ . 

The isoscalar and isovector densities (|3U)) can now be 
expressed by using the auxiliary densities (|4"4"1) - (|4"T)) . and 
they read 

poo(x,x') = cos^(pW(x,x')+p%>\x,x')) 

- ffa£(j>W(x,x')-pW(x,x'j), (48) 
p w (x,x') = cos^(pM(x,x')-pW(x,x')) 

- sm^(p^(x,x') + p^(x,x')), 

pu(x,x') = - V2 cos -p^(x,x') 

+ v/2 sin|pW(x,x') ) 

,5i_i(x,x')= \/2 cos-ffl(x,x') 

+ V2 sin|pW(x,x'). 



D. Hamiltonian kernel 



Conventional Skyrme EDF can be expressed by bi- 
linear forms of six isoscalar (t = 0) and six isovector 
(t = 1) local densities, including the particle p t , kinetic 
r T , spin s t , spin-kinetic T t , current jt, and spin-current 
Jt densities and their derivatives. Standard definitions 
of these densities can be found in numerous references, 
see, e.g., Refs. d, |46[ and references quoted therein. It is 
to be noted that in the standard MF theory the proton 
and neutron s.p. wave functions are not mixed, i.e., the 
proton- neutron symmetry is strictly conserved [46[ [43]. 
Therefore, the MF isoscalar (isovector) densities are sim- 
ply sums (differences) of neutron and proton densities 
since only the third component of the isovector density 
is nonzero. 



The isospin projection essentially preserves the func- 
tional form of the Skyrme EDF derived by averaging the 
isospin-invariant Skyrme interaction over the Slater de- 
terminant. All what needs to be done is a replacement 
of the local density matrix by the corresponding transi- 
tion density matrix. Moreover, the bilinear terms that 
depend of the isovector densities must be replaced by 
the full isoscalar products of the corresponding isovector 
transition densities. Special care should be taken of the 
density-dependent term of the Skyrme interaction as the 
extension of this term to the transition case is undefined 
in the process of averaging the Skyrme interaction over 
the Slater determinant. Following the argumentation of 
Refs. [HI [4§] , we replace the isoscalar density by the tran- 
sition isoscalar density matrix, that is, p n + p p — ► poo- 

As for the Coulomb-interaction kernel, it depends on 
the isoscalar and isovector transition densities in the fol- 
lowing way: 



(49) 



2 A 



e f dri f dr 1 £ £ V^(r u r 2 ), (52) 



A=0 //— — A 



(50) 



(51) where 
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v l±l 



V 2±l 



2±2 



~ |poo(r 2 ,ri)poo(ri,r 2 ) + ~pi(r 2 ,ri) o pi(ri,r 2 ) 
+ soo(T2, r i) • soo(ri,r 2 ) + -Si(r 2 ,ri) o -Si(ri,r 2 ) 



-poo(^i)/5io(r 2 ) + - {poo(r 2 ,r 1 )p w (r 1 ,r 2 ) + s 00 (r 2 , r x ) ■ s w {r 1 ,r 2 )} , 
/5oo(7-i)Pi=fi(»" 2 ) - - {poo(r2,r 1 )p lTl (r ll r 2 ) + s o(r- 2 ,7-i) • s lTl (r 1 ,r 2 )} ■ 



g {pio(n)pio(T-2) + /5n(Vi)/5i-i(r 2 )} 

^ {/5i (r- 2 , ri)pi (ri , r 2 ) + /5n(r 2 , ri)pi_i(ri , r 2 ) 

+ si (r 2) T-i) • §10(7-1, r 2 ) + sn(T- 2 ,n) • si-i(ri,r 2 )} , 

-^Pio(n)pi=Fi(r 2 ) + ^^{pio(r 2 ,ri)pi =F i(ri,r 2 ) + Sio(r 2 ,ri) • si T i(n, r 2 )} , 
^Pi=Fi(ri)pi T i(r 2 ) - ^={^1(^,7-1)^1^1^2) + si =F i(r 2 ,T-i) • §1^1(7*1, r 2 )} . 



(53) 
(54) 
(55) 



(56) 
(57) 

(58) 



Here, the symbols • and o stand for the scalar products 
of vectors and isovectors, respectively. 



E. Coulomb rediagonalization and the isospin 
mixing 

To properly account for the isospin mixing effects, fol- 
lowing Refs. [H, the total Hamiltonian Hnn © 
(strong interaction plus the Coulomb interaction) is redi- 
agonalized in the space spanned by the good-isospin wave 
functions {J), and the resulting eigenstates are denoted 
by 



\n,T z 



T>\T Z \ 



\T,T Z ) 



(59) 



and numbered by index n. The amplitudes T define 
the degree of isospin mixing. In particular, the isospin- 
mixing parameter for the lowest energy solution -E„ = i ; t 2 
is defined as ac = 1 — la^j 1 ^! 2 . 

Precise determination of the Coulomb mixing con- 
stitutes a notoriously difficult problem, but is strongly 
motivated by its impact on fundamental physics tested 
through the super-allowed /? decays [5(1 HH, H3] • Of im- 
portance here is to capture the proper balance between 
the short-range strong interaction and the long-range 
Coulomb polarization. This balance is naturally taken 
into account by the DFT approach, but is not accessi- 
ble within a perturbative-analysis theory [53| or hydro- 
dynamical model [HI . This is illustrated in Fig. [3] that 
shows the excitation energies of the doorway T=l states 
in e-e N=Z nuclei, i.e., the energies of the n = 2 states 
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FIG. 3: Excitation energies of the doorway T=l states in the 
e-e N=Z nuclei relative to the SHF g.s. energies Erf- Dia- 
monds, dots, and circles show the self-consistent results ob- 
tained in the isospin-projected DFT with SIII [4l]], SLy4 [551 ]. 
and SkP [56| Skyrme functionals, respectively. Horizon- 
tal lines mark the mean DFT values. Thick dotted lines 
show estimates based on the perturbation theory, Et=i ~ 
82/A 1/3 MeV [H [HI], and on the hydrodynamical model, 
E T =i » 169/A 1/3 MeV 0. 



in Eq. (]59[) relative to the HF g.s. energies Ehf- The 
estimates given by the perturbative/hydrodynamical ap- 
proaches are clearly at variance with the self-consistent 
results. 

Figure [3] shows that not only the values but also the 
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FIG. 4: Superdeformed bands in Ni. Experimental data 
for Band 1 (dots) and Band 2 (circles) of Ref. [20I ] are com- 
pared with HF+SLy4L predictions marked by full and open 
triangles, respectively. The left panel shows the standard HF 
results while the predicted T=0 SD bands obtained in the 
isospin-projected approach are displayed in the right panel. 



A-dependence of the doorway excitation energies differ 
substantially from the self-consistent results, pointing 
very clearly to a non-perturbative origin of the Coulomb 
mixing. One can also notice that the mean excitation 
energies of the doorway states (directly impacting the 
Coulomb mixing) strongly depend on the EDF param- 
eterization. The currently used Skyrme functionals are 
not sufficiently constrained in the isospin sector to pro- 
vide reliable estimates of ac ■ At this point is also not at 
all obvious what parts of EDF have to be refined in or- 
der to improve on this situation, see Ref. [29[ for further 
details. 



IV. ILLUSTRATIVE EXAMPLES 

All calculations presented in this work were obtained 
by using isospi n p ro j ection implemented in the DFT 
solver HFODD [57|,[58|. We used the harmonic-oscillator 
basis of iVn=12 shells. We considered SIII {41], SLy4 H, 
and SkP [56[ Skyrme functionals, as well as for their Lan- 
dau and tensor versions denoted by subscripts L and T, 
respectively, see Ref. [59| for details. Pairing was ne- 
glected. 



A. Superdeformed bands in Ni 

We begin discussion from the case of SD bands in 56 Ni. 
As discussed in Sec. HU the standard cranked HF theory 
gives a reasonable reproduction of Band 1 while there are 
problems with theoretical interpretation of Band 2 built 
on p-h excitations to the [440] 1/2 proton and neutron 
levels. The results are summarized in Fig. 2) The HF 
(left panel) and isospin-projected (right panel) calcula- 
tions were performed with SLy4/, functional. Since the 
HF configuration corresponding to Band 1 is predomi- 
nantly isospin symmetric, it is only weakly affected by 
the isospin projection. The calculated isospin impurity 
within this band is small (~ 2%). This value is similar to 



that obtained for the spherical g.s. configuration of 56 Ni. 
On the other hand, the Band 2 based on the 7r[440]l/2 
orbital represents, before rediagonalization, almost equal 
mixture of the T—0 and T—l components. Therefore, its 
isospin impurity assumes an unrealistic value of about 
50%. 

Following the isospin projection, the low-spin part of 
the T=0 SD band originating from Band 2 is shifted down 
by about 0.95 MeV with respect to the unprojected HF 
result. The isospin impurity in this band is in the range 
of 6% to 8%, and it slowly increases as a function of the 
angular momentum. In fact, these values may indicate a 
presence of an uncontrolled isospin mixing related, most 
likely, to the angular-momentum non-conservation. As 
shown schematically in Fig. [2j a T—l SD partner struc- 
ture is expected to lie higher in energy. We indeed cal- 
culate the T=l band (not shown in Fig. 2]) to lie about 
2.5 MeV above the T=0 band. As expected, the T=0 
and T=l bands projected from the j/[440]1/2 MF config- 
uration are almost identical to those projected from the 
7r[440]l/2 MF configuration. 

By comparing the results of isospin-projected calcula- 
tions with experiment, we see a significant improvement 
as compared to standard HF. The T=0 Band 2 agrees 
well with experiment, both in terms of excitation energy 
and Mol. The predicted crossing between T=0 Bands 1 
and 2 occurs around spin 14 h, i.e., ~ 2h too high as 
compared to the data. This discrepancy, however, rather 
reflects a deficiency of our model in describing the Mol 
of Band 1, which is slightly overestimated as our calcu- 
lations neglect pairing correlations that are expected to 
be important in this configuration (2p| . The large energy 
splitting between the T=0 and T=l Band 2 doublet is 
consistent with observation of only one SD side band. In 
short, the isospin-projected MF theory provides a quan- 
titative explanation of experimental data on collective 
band structures in 56 Ni. 



B. Isospin symmetry breaking at band termination 

Terminating states or seniority isomers are fully- 
stretched p-h configurations with a maximum-spin I max 
that can be built within a given SM space of valence par- 
ticles. Because of their simple SM character, terminating 
states provide a robust probe of SM and MF theories and 
corresponding effective interactions. In this context, of 
particular interest are the terminating states associated 
with the [/ 7 ™ 2 ]/ mra and [d^, 2 f^ ]j mM configurations (in 
the following, n denotes a number of valence particles 
outside the 40 Ca core) in 20 < Z < N < 24 nuclei from 
the lower— fp shell (^4 ~ 44), which were systematically 
measured during the last decade [6fJ, |6l|, [62j, |63|, [6J, |65j . 

According to MF calculations, these specific states ap- 
pear to have almost spherical shapes; hence, the cor- 
relations resulting from the angular-momentum restora- 
tion are practically negligible there (66j. Hence, they 
can be regarded as extreme cases of an almost undis- 
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FIG. 5: Energy differences AE (|60[) between the terminating 
states in Ca,Sc,Ti,and Cr N—Z nuclei. Open and full trian- 
gles mark, respectively, standard and isospin-projected HF re- 
sults for the [dg/j/^t'lL 



(up triangles) and [^/"/VlL 



(down triangles) configurations. Calculations were performed 
by using the SIII_t functional [SjJ. Experimental data (dots) 
and SM results of Ref. [6j| (circles) are shown for comparison. 
Note the limited energy scale. 



turbed s.p. motion, thus offering an excellent playground 
to study, amon g o thers, time-odd densities and fields, 
spin-orbit force [67], [H, [6!| , tensor interactions [7(| , and 
the isospin dependence of cross-shell (sd— fp) p-h matrix 
elements [68l. l7lj|. 

While most of the terminating states are uniquely de- 
fined, the [rf^/^ li'm^ states in the N=Z nuclei pro- 
vide a notable exception. Indeed, within the MF approxi- 
mation these particular states can be created by promot- 
ing either one proton, [^3/2/7/2]^ ; or one neutron, 
K3/2 ^7/2 1 ] / m0 x' across the magic gap 20, Because the 
Coulomb energy difference between these configurations 
is small, the resulting energy levels are almost degener- 
ate: E^f^Y^J * W/l^lLJ' Therefore, 
we encounter problems related to isospin symmetry sim- 
ilar to those discussed earlier in Sec. [Til 

As in the case of 56 Ni, we encounter the situation where 
the isospin symmetry is manifestly broken by the MF 
approximation and the predictions are at variance with 
empirical data. The difficulty with the isospin content 
of terminating states was recognized in Ref. (f38| . where 
a purely phenomenological method of isospin restoration 
was proposed. It has resulted in a good MF description 
of experimental data, at the level of the state-of-the-art 
SM calculations. Quantitatively, however, the estimated 
energy correction due to the isospin projection appeared 
to have surprisingly strong A-dependence, changing quite 
rapidly from 5E T rz2MeV in 40 Ca down to SE T falMeV 
in 46 V, see Fig. 4 of Ref. 68]. This trend has been found 
to depend weakly on the EDF parameterization. 



In the present work, we repeat calculations of Ref. [68| , 
however, by using the mathematically rigorous isospin 
projection of Sec. IIIIl The energies of the isospin- 
projected states [^3/2/7/2 ^f^ax are shown in Fig. [5] rela- 
tive to those of the [f?/ 2 ]i maa: configurations: 



AE 



7/2J 
^([^3/2/7/2 



,aJ - E ([f7/2\l ma J- 



(60) 



Full up and down triangles show results obtained by pro- 
jecting the T=0 component of the Slater determinant 
corresponding to either one proton, I [^3/2 /7/2 1 ] ) ' or 
one neutron, | [^3/2/7/2 ]i maw )i P" n excitation through the 
magic gap 20, respectively. For comparison, the unpro- 
jected HF energies for the proton and neutron configu- 
rations are shown with the open up and down triangles, 
respectively. Here, all calculations were performed by 
using the SIII [4l[ Skyrme functional. 

The isospin-projected results obtained from the 

fn+lli/ 



IL U 3/2J7/2 Jl„ 



and |[^3 / / 1 2/7/2 i ]/„„ T ) HF configurations 



are similar but not identical, reflecting small polarization 
differences due to proton and neutron p-h excitations. 
Note, that slightly better results are obtained by pro- 
jecting from the proton configurations as the resulting 
levels appear slightly lower in energy. This is probably 
not surprising, because these states include directly the 
polarization of the Coulomb field by the proton p-h exci- 
tation and thus they have slightly richer isospin structure 
than their neutron counterparts. 

The results of rigorous isospin projection closely fol- 
low those obtained in the phenomenological approach of 
Ref. (68[. As shown in Fig.O the energy corrections cal- 
culated in the isospin-projected HF, 



5E 7 



E{[d, 



— 1 £Ti 
3/2J7/2 



n+i-|T=l 

/ma* 



^([^3/2/7/2 J/, 



liT=0 



(61) 

also exhibit an appreciable A-dependence, in a close anal- 
ogy to the previously obtained phenomenological trend. 
Consequently, the A-dependence of AE ([BTJ)) , calculated 
using the isospin projected EDF approach also shows a 
different pattern than experiment and SM. This result 
is weakly dependent on EDF parameterization, suggest- 
ing that there is some generic problem pertaining to the 
standard form of the Skyrme EDF. 

Recently, we have developed a new class of the Skyrme 
functionals with spin-orbit and tensor terms locally refit- 
ted, to reproduce the /V/2 — /s/2 spin-orbit splitting in 
40 Ca, 48 Ca, and 56 Ni [§J. The spin-orbit and the tensor 
strengths obtained in this way turned out to be fairly 
independent of other coupling constants of the Skyrme 
functional. This result indicates that the strong depen- 
dence of the spin-orbit strength on EDF parameteriza- 
tion, in particular on the isoscalar effective mass [z3], is 
likely to be an artifact of fitting protocols based predom- 
inantly on data pertaining to bulk nuclear properties. 
Indeed, as discussed in Ref. |73| , the use of inaccurate 
models in the fitting procedure can lead to results that 
strongly depend on the fitting protocol itself; hence, can 
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FIG. 6: Full symbols label energy correction SEt (|61[) cal- 
culated by projecting good isospin from the HF states corre- 
sponding to the [^ 1 2 / 7 "^ 1 ]/ mra configurations in N=Z nuclei 
and subsequent Coulomb rediagonalization. Open symbols 
represent results of a phenomenological method proposed in 
Ref. [(3|| . Calculations were carried out by using the SIIIl 
(dots) and SkP.LT (triangles) Skyrme functionals [59| . See 
text for further details. 



result in contradictory information on the key model pa- 
rameters (compare, e.g., results of Refs. [72[ and [z3|)- ^ 
seems that this is exactly the case for the current param- 
eterizations of the Skyrme EDF. As shown recently in 
Ref. [75(, parameterizations that correctly describe the 
spin-orbit properties in light nuclei do not fare well in 
heavier systems. This again points to limitations of the 
second-order Skyrme EDF [76[ and to a danger of draw- 
ing^ conclusions on tensor interactions from global fits 

Applications of new functionals to the terminating 
states in N ^ Z, A ~ 44 nuclei [ztJ have revealed that 
removal of the artificial isoscalar effective mass scaling 
from the spin-orbit restores the effective mass scaling in 
the s.p. level density. As a consequence, only the forces 
having large isoscalar effective masses {Mj- > 0.9) such 
as SkPr and SkOr [2Q[, are able to reproduce empirical 
data involving s.p. levels in light nuclei. 

This conclusion is nicely corroborated by results pre- 
sented in Fig. [3 which shows the predictions for ter- 
minating states using the new Skyrme parameteriza- 
tion SkPiT- R is rewarding to see that the modified 
functional yields results that are very consistent with 
SM. This is particularly true for the isospin projection 
from the | [^3/2/7/2 ]? mM ) configurations. The new val- 
ues of AE reasonably agree with experiment, consid- 
ering the energy range of the plot, While the detailed 
A-dependence is reproduced very well, theoretical AE 
curve slightly underestimates experiment. Identification 
of a specific source of this remaining discrepancy requires 
further studies. 
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V. SUMMARY AND CONCLUSIONS 

In this work, we carried out the theoretical analysis 
of isospin breaking in nuclei around N=Z based on the 
density functional theory. We show that the sponta- 
neous breaking of isospin symmetry inherent to MF lim- 
its the applicability of self-consistent theories, such as HF 
and DFT, to nuclear states with T—0. To remedy this 
problem, we propose a new isospin-symmetry restoration 
scheme based on the rediagonalization technique in good- 
isospin basis. 

The isospin projection algorithm is described in detail, 
including the derivation of essential Hamiltonian kernels. 
We demonstrate that isospin projection is free from prob- 
lems plaguing other symmetry restoration schemes (par- 
ticle number, angular momentum) related to kernel sin- 
gularities. 

Applications of the isospin-projected DFT approach 
include the SD bands in 56 Ni and terminating states in 
N=Z medium-mass nuclei. In both cases, the symmetry 
restoration remedies previously noted deficiencies of the 
HF method and significantly improves agreement with 
experiment. The results for terminating states corrob- 
orate previous suggestions that large isoscalar effective 
masses are needed to reproduce experimental spectro- 
scopic data involving s.p. levels. 

The examples presented in this work primarily concern 
high-spin configurations in medium mass nuclei in which 
pairing correlations are expected to be weak. To be able 
to address theoretically a variety of low-spin phenomena 
and observables around the N=Z line, such as ground 
states of odd-odd nuclei, binding energy staggering, su- 
perallowed beta decays, and charge-exchange reactions, 
isovector and isoscalar pairing correlations must be incor- 
porated and the proton-neutron symmetry of HF must be 
broken before isospin projection j4(|[47|]. Those will be 
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subjects of our forthcoming studies. 
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